Sheared active fluids: thickening, thinning and vanishing viscosity 
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We analyze the behavior of a suspension of active polar particles under shear. In the absence of 
external forces, orientationally ordered active particles are known to exhibit a transition to a state of 
non-uniform polarization and spontaneous flow. Such a transition results from the interplay between 
elastic stresses, due to the liquid crystallinity of the suspension, and internal active stresses. In the 
presence of an external shear we find an extremely rich variety of phenomena, including an effective 
reduction (increase) in the apparent viscosity depending on the nature of the active stresses and the 
flow-alignment property of the particles, as well as more exotic behaviors such as a non-monotonic 
stress/strain-rate relation and yield stress for large activities. 



I. INTRODUCTION 

Colonies of swimming bacteria, in vitro mixtures of cy- 
toskeletal filaments and motor proteins, and vibrated 
granular rods are examples of active systems composed 
of interacting units that consume energy and collectively 
generate motion and mechanical stresses. Due to their 
elongated shape, active particles can exhibit oricntational 
order at high concentration and have been likened to "liv- 
ing liquid crystals" [T]. Their rich collective behavior in- 
cludes nonequilibrium phase transition and pattern for- 
mation on mesoscopic scales. It has been modeled by con- 
tinuum equations built by modifying the hydrodynamics 
of liquid crystals to include nonequilibrium terms that 
account for the activity of the system [2H1]j or derived 
from specific microscopic models [3 [6] . 

A striking property of confined active liquid crystals is 
the instability of the uniform aligned homogeneous state 
and the onset of spontaneously flowing states, both sta- 
tionary and oscillatory |8] . This occurs because local 
orientational order generates active stresses that are in 
turn balanced by flow, yielding a state that can support 
local inhomogeneities in the flow velocity and the local 
alignment, while maintaining a net zero force. Loosely 
speaking, a confined active liquid crystal "shears itself" 
even in the absence of externally applied forces. It is 
then not surprising that the rheology of such active liq- 
uid crystals in response to an external shear will be very 
rich. 

Phenomenological work by Hatwalne and collabora- 
tors [9] first pointed out that activity lowers the linear 
bulk viscosity of tensile suspensions, such as most swim- 
ming bacteria, while it enhances the viscosity of con- 
tractile systems, and that this enhancement may become 
very large near the isotropic-nematic transition. A semi- 
microscopic model of contractile suspensions of motor- 
filaments mixtures confirmed these results and predicted 
an actual divergence of the viscosity of contractile sus- 
pensions at the transition [10] . Recent numerical studies 
of active nematic films by Cates et al. |llj have confirmed 



that this result survives when the effect of boundaries is 
included. In addition, it was found that tensile nematic 
suspensions can enter a regime of vanishing apparent vis- 
cosity in proximity of the isotropic-nematic phase tran- 
sition. Such a "superfluid" window was interpreted by 
the authors of Ref. [TT] as the appearance of bulk shear 
bands accommodating a range of macroscopic shear-rates 
at zero stress. Finally, the predicted activity-induced 
thinning of bacterial suspensions has been demonstrated 
in recent experiments in Bacillus subtilis |12H14) . 

Active particles exert forces on the surrounding fluid, 
resulting in local tensile or contractile stresses propor- 
tional to the amount of orientational order, afj ~ anirij , 
where a is proportional to the force exerted by the active 
particles on the fluid and n a unit vector denoting the di- 
rection of broken orientational symmetry. The sign of a 
determines whether the flow generated by the active par- 
ticles is tensile (a < 0) or contractile (a > 0). In the case 
of swimming organisms, the former situation describes 
"pushers", i.e., most bacteria (e.g., E. Coli), while the 
latter corresponds to "pullers" (e.g., Chlamydomonas) 
(see Fig. [T]). An important distinction between uniaxial 
active particles concerns the possibility of forming phases 
with or without a non-zero macroscopic polarization. Ap- 
olar particles are fore-aft symmetric and can form ne- 
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FIG. 1: (color online) Schematic example of the flow field 
surrounding a tensile (left) and contractile (right) swimming 
microorganism. 
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matic phases in which macroscopic quantities are invari- 
ant for n — > n. Polar particles can also form phases 
characterized by a non-zero macroscopic polarization in 
the direction of a polar director p in which they undergo 
collective motion with mean velocity v ~ (3 p, with /3 
is the typical self-propulsion velocity. This directed mo- 
tion occurring in polar suspensions contributes to a non- 
equilibrium local stress of the form cr^ ~ j3 (dipj + djPi). 

Most theoretical work has focused on the rheology of 
active nematic (j3 = 0), while the shear response of ac- 
tive polar suspensions is far less explored [HI [T3] . We 
find that for a fixed value of /3, the behavior of active 
suspensions depends on the interplay between the local 
contractile/tensile stresses, embodied in the parameter a, 
and the flow-aligning behavior of liquid crystalline par- 
ticles, described by the flow alignment parameter, A [? 
]. Rod-shaped particles typically have A > 0, spherical 
particles have A = 0, while the case A < describes disk- 
shaped molecules such as those found in discotic liquid 
crystals. In passive liquid crystals the magnitude of A 
controls how the director field responds to a large shear 
flow away from boundaries. For |A| > 1 the director tends 
to align to the flow direction at an angle 6q such that 
cos 2^o = 1/A, while for |A| < 1 it forms rolls throughout 
the systems. These regimes are known as "flow-aligning" 
and "flow-tumbling" respectively. Understanding of the 
complex rheology of polar and nematic active suspensions 
requires exploring the full parameter space, including the 
important role of boundary conditions. One of the im- 
portant results of this work is a remarkable exact duality 
that holds in the regime where the stress-strain relation 
is linear and shows that tensile (a < 0) rod-shaped flow- 
aligning particles (A > 1) are rheologically equivalent 
to contractile (a > 0) discotic flow-tumbling particles 
( — 1 < A < 0). Using this result, we present below a 
unified description of the linear rheology of active sus- 
pensions of both polar and apolar particles. Some of the 
results are summarized in the "phase diagram" of Fig. [2] 
This figure shows that the rheological properties of an ac- 
tive film subject to an external shear are closely related 
to the onset of spontaneous flow in the absence of shear, 
highlighting the parallel role played in active system by 
mechanical driving forces, such as a macroscopic strain 
rate, and internal active driving forces proportional to a 
and (3. 

An unsheared active film exhibits a transition from the 
homogeneous aligned state to a "spontaneously flowing" 
state, characterized by spatially inhomogencous velocity 
and director profiles [7] . The transition occurs at a criti- 
cal activity a c \ in a film bounded by one no-slip substrate 
and a surface that can freely slide, and at a larger value, 
«c2 > oc c \, in a film bounded by two no-slip planes. The 
lines separating regions of different shades in Fig. [2] are 
the boundaries a c i(/3, A) [see Eq. Q below] separating re- 
gions of spontaneous flow (|a| > a c i) from regions where 
the homogeneous aligned state is stable (\a\ < a c \). In- 
terestingly, when the film is subject to an external shear, 
we find that the flow properties change their qualitative 
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FIG. 2: (color online) The figure displays the regions of pa- 
rameters where spontaneous flow occurs in an unsheared ac- 
tive film on a substrate. The regions of spontaneous flow 
are bounded by the critical activity a c i(/3) given in Eq. Q 
(solid and dashed lines) and are shaded orange, with lighter 
shades corresponding to increasing values of j3. The same 
critical activity also separates the regions |a| < a c i where 
the theoretical stress-strain curves are monotonic and the ac- 
tive suspension is either thinned or thickened by activity at 
small shear rates, as indicated, from the regions \a\ > ot c \ 
where the theoretical stress-strain curves are nonmonotonic, 
with possible "superfluid" or hysteretic behavior. 



behaviour at exactly these same critical values of activ- 
ity. For a c \ < \a\ < a C 2, the theoretical stress-strain 
rate curves obtained from our one dimensional model are 
nonmonotonic (see Fig. [8]) and the active suspension is 
strongly non-Newtonian. We suggest a number of dif- 
ferent interpretations of the nonmonotonic part of the 
stress-strain rate curve shown in Fig. [7 These include 
macroscopic "superfluid-like" behaviour TT] with zero ef- 
fective viscosity, yield-stress behaviour or hysteresis. Fi- 
nally, for \a\ > a C 2, the theoretical stress-strain curve has 
a discontinuous jump at zero strain rate, corresponding 
to a finite "spontaneous stress" in the absence of applied 
shear [10] . 



II. THE MODEL 

Our model of active suspension consists of a two- 
dimensional film of rod-like particles of length £ confined 
to a channel of infinite length along the x axis and finite 
thickness L along the y axis (see Fig. [3]). Because of 
the chosen geometry, the system is invariant for trans- 
lations along the x axis. The total density of the sus- 
pension, p = Mc + Psoivcnt, with c the concentration of 
active particles and M their mass, is assumed to be con- 
stant, thus V • v = 0, with v the flow velocity. We as- 
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FIG. 3: (color online) Schematic representation of a quasi- 
one-dimensional film of thickness L. In our model the film is 
sitting on a non-slipping surface and is sheared from the top 
at constant velocity vo- The polar rods form an angle 9 with 
respect to the infinite direction x of the film. Because of the 
quasi-one-dimensional geometry, the system is invariant for 
translations along the x axis. 



sume that the film is sheared at a constant (macroscopic) 
rate 7 by keeping the lower plate at y = fixed, while 
the upper plate at y = L is moved at constant veloc- 
ity vq. The macroscopic shear-rate is defined then as 
7 = Vq/L = J Q (dy/L)u, where the rate-of-strain ten- 
sor Uij = (diVj + djVi)/2 has only non-zero components 
u xy — u yx — d y v x /2 = u/2. Theoretical stress-strain 
curves are obtained by fixing the macroscopic strain rate 
7 and calculating the resulting stress a. 

We consider a polarized active suspension and focus 
only on spatial variations in the direction of the polar- 
ization P. The hydrodynamic equations for an active 
polar suspension have been formulated by incorporating 
the active contributions (proportional to the rate of en- 
ergy consumed by the active units) into the hydrody- 
namic equations of a passive polar liquid crystalline film. 
Some of the active contributions, discussed above, are 
not allowed by the conditions which define liquid crystal 
systems at equilibrium and hence are intrinsic to active 
systems. Other terms have the same form as those of 
passive polar liquid crystals and can simply be included 
by modifying the prefactors of the terms obtained from 
a passive systems. As such, the modified "passive" con- 
tributions to the equations of motion can be described 
starting from the non-equilibrium analogue of the Frank 
free-energy of a suspension of polar particles in a solvent: 
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with C the compressional modulus and K\ and K3 
the splay and bend elastic constant. The parameters 
di,Bi,Ki,C are understood to have both passive and ac- 
tive contributions. In the following we will take K\ = 
K3 = K. The last three terms in the expression of the 
free-energy couple concentration and splay and are also 
present in equilibrium polar suspensions. 



The dynamics of the concentration and the polariza- 
tion are described by 



d t c 



-v • [ c (v + cftP)+r'h + r"f] 



(la) 



[8 t + {v + cP3P)'V\P i + u iJ P J 

= Xu ij P j +Th i +T'f i , (lb) 

with Uij — (diVj — djVi)/2 the vorticity tensor, h = 
-8F/SP the molecular field and f = -V{5F/8c). The 
flow velocity satisfies the Navier-Stokes equation [T5] : 



p(d t + v • V)vi 



(2) 



with V • v = to guarantee incompressibility, and stress 
tensor given by dissipative, reversible and active contri- 
butions, <Tij = 2r\u v} + a r i} + &fj + , with 
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[diPj + djPi + SijV'P] 
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(3b) 



o£ - -1% - -(Pihj + P s hi) + ~(Pihj - Pjhi) , (3c) 

where II is the pressure, 77 the shear viscosity, and we 
have assumed an isotropic viscosity tensor. We now con- 
sider a solution deep in the polarized state and neglect 
fluctuations in the magnitude of the polarization, i.e., 
assume |P| = \J —aij a^. For simplicity we also redefine 
units so that |P| = 1. The condition P — constant deter- 
mines the longitudinal part h\\ = p ■ h of the molecular 
field that can then be eliminated from the hydrodynamic 
equations. The details associated with imposing the con- 
stancy of the magnitude of the polarization and deriving 
the hydrodynamic equations solely in terms of the polar 
director p = P/|P| are given in Appendix |A} With this 
choice, the hydrodynamic equations for p and c can be 
written in the form 

9 t c+ V-c(v + /?!cp) = di [DijdjC + Xj'ukipkpip,} , (4a) 



[d t + (v + /3 2 cp) • V]p 4 + Uijpj 



^UjkPk H die (3,-V • p + n\7 2 pj 

c c 



, (4b) 



with 7' = ryr, « = tk, w = t(b 1 - b 3 ) and s? 3 = 

Sij — PiPj the transverse projection operator. is an 
effective diffusion tensor given by 

Dij = DxSij + D 2 piPj , (5) 

where D\ = D — 7'w/co and D2 = "/'w/cq — D£. Finally, 
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TABLE I: Estimates of active parameters for two types of 
active suspensions: (i) mixtures of cytoskeletal filaments and 
cross-linking motor proteins [5j |6j [16] , with m a dimensionless 
density of crosslinking motor clusters , Mq the speed at which 
motor proteins walk on filaments, in turn proportional to the 



rate of ATP consumption, and 



u t 



with £ m the size 



of a motor cluster; and (ii) swimming microorganisms [15] , 
where v sp ~ (//C) 6 is the self-propulsion speed of an individ- 
ual organisms, with / the force that swimmers exert on the 
fluid, e < 1 a dimensionless number determined by the shape 
of the swimmer and C, ~ 1/F is the longitudinal friction coeffi- 
cient of a rod-like swimmer of length £. For both systems the 
precise values of parameters obtained from each microscopic 
model differ from the above by numericsl constants of order 
unity. 



the reversible part of the stress tensor of - becomes: 
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The equations for an active suspension have been written 
down phenomenologically and also derived from various 
semi-microscopic models. The structure of of the equa- 
tions is generic and applies to a broad class of "living 
liquid crystals" . The parameters in the equations are 
of course system and model specific. In motor/filament 
mixtures activity arises from clusters of motor proteins 
crosslinking pairs of filaments. The active couplings are 
therefore of order c 2 in this case [6] . In suspensions of 
swimming microorganisms, activity can be described in 
terms of the active force / that each swimmer exerts on 
the surrounding fluid. In this case the active couplings 
arise even at the single-swimmer level and are of order 
c [15 . Estimates for the active parameters obtained from 
semimicroscopic models are summarized in Table [TJ The 
equations for an active nematic can be obtained from 
those of a polar systems by setting /?, = w = 0. In the 
following we assume /3i = /3 3 = — /3 2 — ft, as appropriate 
for motor filament-systems. 

It is convenient to work with dimensionless quantities. 
Spatial variables are normalized with the length I of the 
rods. Thus y — >■ y/i. Temporal variables are normalized 



with the time scale of splay and bending fluctuations, 
thus t — > tjr where r — l 2 /n. A mass scale is set by 
r/r. All the other quantities are normalized accordingly. 
In these units the hydrodynamic equations for the rods 
concentration (f> = c/cq, with Cq the mean density, and 
the director/polarization angle 9, with p = (cos 9, sin 9), 
for the geometry of interest are 



p(d t + Vydy)V x = dy(J X y 



(6a) 



dtcb = d. 



sm0 + V(0)d y (j) + Au sin (2 sin 2(9} , (6b) 



d t 9 = -/ 



i cos 9d, 



dp + jc{e)d 2 y t 



+wcos6sin9(d y 9) 2 - u(l - Acos26>) , 



(6c) 



where V{9) = £>(!-£ sin 2 



; cos 2 9 is a diffusion co- 



efficient, JC(9) = 1 — wcos 2 9 describes the energy cost of 
bend and splay deformations, and A is the flow- alignment 
parameter. In a steady state the stress tensor <j xy = a is 
constant across the film and it is given by 



■q + A 2 sin 2 29 + Xw sin 2 9 sin 29 (d y 



+ [w — Xw — X(w — wo) cos 29] cos 9d y <f) 



■a(j) 2 sin20-2^> 2 sin 9d v 9 , 



(7) 



with r\ the bare viscosity and Wq a constant proportional 
to the ratio between the translational and orientational 
diffusion coefficients (i.e. wq ~ D/K). Our goal is to 
study the relation between the induced shear stress a 
and the applied shear rate 7 as a function of the two 
fundamental active parameters a and /3 representing the 
magnitude of the internal contractile/tensile stress and 
the velocity scale of directed motion. In order to con- 
struct a a vs 7 map, we integrate Eqs. Q numerically 
with boundary conditions v x (0) = and v x {L) = vq, 
9(0) = 9(L) = and j y (0) = j v (L) = which im- 
plies 4>'(0) — 4>'{L) — 0. As initial conditions we choose 
%,0) = 0and <My,0) = l. 

In the absence of applied shear, active polar and ne- 
matic films exhibit a transition from a quiescent (v x = 0) 
aligned (9 = 0) state to a state of spontaneous flow, with 
both inhomogeneous alignment and velocity profiles. The 
critical value of activity where the instability occurs de- 
pends on boundary conditions. For a film bounded by a 
no-slip substrate and a surface that can freely slide it is 
given by [5]: 



a cl (/3, X) 



r)(l-w) Pw[rj + (1 - A) 2 ] 
2$(1 - A) + 2(1- X){D-w) 



(8) 



and the spontaneously flowing state has a — 0. For a 
film bounded by two no-slip surfaces the critical value is 
a c2 = 4a cl and the spontaneously flowing state is charac- 
terized by a finite value of a. The regions of spontaneous 
flow in the (A, a) plane are displayed in shades of orange 
111 Fig. [2| In these regions the film exhibits strongly non- 
linear rheology, with nonmonotonic stress-strain curves, 
as described below. 
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FIG. 4: Stress (er) vs strain (7) for an active nematic (J3 = 
w = 0) suspension for various a. Flow-tumbling system with 
A = 0.1 are marked by circles and flow-aligning systems with 
A = 1.9 by triangles. Other parameters are set L/£ — 5, 
rj = 1, <j>o — 1, D — 1 and £ = 0.3. The inset shows the 
comparison with the analytical result given in Eq. |9J. 



III. LINEAR RHEOLOGY OF WEAKLY 
ACTIVE SYSTEMS 



For |a| < a c i, corresponding to the gray regions of Fig. [2] 
the stress strain curves are monotonic and remain lin- 
ear over a broad range of 7, as shown in Fig. [4] Non- 
Newtonian behavior sets in at smaller values of 7 with 
increasing a. As the value of a is increased the slope of 
the linear portion of the stress-strain curves for a < a c \ 
decreases with increasing a, indicating that contractile 
active stresses lower the effective viscosity of the system. 
The effective linear viscosity can be calculated analyti- 
cally by solving Eqs. (6c) and (6b) perturbatively in a 
by expanding the fields 9 and (f) &s 9 = 9 + o9\ + a 2 9 2 • • • 
and </> = (f>o + a(f)i + a 2 9 2 ■ ■ ■ The quantities #0 and </>o rep- 
resents here the stationary solution of the hydrodynamic 
equations in absence of shear flow. If the suspension is in 
an aligned state at t = 0, when the shear is switched on, 
then 9q — and <fio = const. We note, however, that this 
perturbation analysis breaks down in the region a > a c \ 
of spontaneous flow, as in that case both 6, <fi are spatially 
varying even at a — 0. It is straightforward to solve Eqs. 
(6c) and (6b) to first order in a. We then obtain the lin- 
ear apparent viscosity defined as 77 app = lim^o er/7 and 



given by 



where tanc(ir) = tan(x)/ir. and 
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(10a) 
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For passive system a — (3 = w = 0, and 77 app = ry, 
as expected. For active nematic, /3 = w = and the 
apparent viscosity is simply 



tanc 



k L\ 
2 11 



(11) 



with k — ^/2a</>o(l — X)/i]- K a (l — A) < 0, fc is imag- 
inary and the tan function at the denominator of ?7a PP 
is replaced by its hyperbolic counterpart. Since tanh(x) 
increases more slowly than x, the resulting apparent vis- 
cosity will increase. If a(l — A) > 0, k is real and since 
the tan(a:) function grows more rapidly than x we ex- 
pect then a rapid decrease in the apparent viscosity as 
I a I is increased. This shows that the linear rheology of 
pullers/contractile systems with A < 1 are the same as 
those of pushers/tensile systems with A > 1. From Eq. 
([9]) it is indeed simple to prove that the apparent viscosity 
?7app is invariant under the transformation 



f?a PP (a, /?, A) = 77a PP (-a, /3, 2 - A) . 



(12) 



Thus flow-aligning pullers with A = 1 + e (for < 
e < 1) will exhibit the same apparent viscosity of ow- 
tumbling pushers with A = 1 — e: 7y app (— |a|, (3, 1 + e) = 
?7app(M, ft, 1 — e). This duality is displayed in the top 
frame of Fig. [6] that shows the linear apparent viscosity of 
active nematic suspensions as a function of |cv| for several 
values of A. The solid curves (red online) show that both 
contractile/flow tumbling suspensions and tensile/flow 
aligning ones are thinned by activity. The dashed curves 
(blue online) refer to either contractile/flow aligning sus- 
pensions or tensile/flow tumbling ones and show that 
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FIG. 5: Schematic example of the flow field surrounding a 
tensile/flow-aligning (right) and contractile/flow-tumbling ac- 
tive particle. For the choice of the parameters a and A given 
in Eq. ( 12 I the two flows are identical, leading to an equal 



apparent viscosity. 



G 







'J 

J-.-9 ; 


1.5 












1 












0.5 
















\0.1 


0.0 







-0.2 -0.1 0.0 0.1 0.2 







0.1 


"{.85 






0.3 








0.5 


1.65 













2 


3 


4 





FIG. 6: (color online) Apparent viscosity 77 app for active ne- 
matic (top) and polar (bottom) suspensions. Solid/red lines 
represent flow-tumbling systems (A < I) while dashed/blue 
lines flow-aligning systems (A > 1). The corresponding val- 
ues of A are indicated next to the lines. In the bottom plot 
a was set to zero. The top frame emphasizes the duality 
discussed in the text. 



these systems are thickened by activity. Bacteria such 
as E-Coli are pushers (a < 0) and generally elongated 
in shape, corresponding to A > 1. Our results therefore 
confirm the activity-induced thinning of bacterial suspen- 
sions first predicted by Hatwalne et al [9] and recently 
observed in [13] . In contrast, algae like Chamydomonas 
that propel themselves from the front (and are therefore 
pullers, with a > 0). Whether they are thickened or 
thinned by activity depends intimately on their shape, 
i.e. on whether they can be described as objects with 
A > 1 or A < 1. Similarly, motor/filament mixtures 
are generally contractile (a > 0) are are expected to be 
thickened or thinned by activity depending on the effec- 
tive value of A. 

This duality has a simple interpretation. Active con- 
tractile (tensile) particles produce an ingoing (outgoing) 
flow in the surrounding fluid, but while flow-aligning par- 
ticles orient at a positive angle with respect to the flow 
direction, flow-tumbling particles orient at a negative an- 
gle under a small applied shear (see Fig. [5j). As a re- 
sult, the average flow fields produced in the surround- 
ing fluid are identical in the two cases and produce the 
same resistance to the imposed shear flow. This equiv- 
alence holds only for small applied shear stresses. For 



large shear-rates the configuration of the director field of 
a flow-tumbling suspension is dramatically different from 
that of flow-aligning one and the similarity between the 
two flow- fields no longer holds. 



IV. NONLINEAR RHEOLOGY OF STRONGLY 
ACTIVE SYSTEMS 

The linear apparent viscosity given by Eq. ( |lTj ) van- 
ishes at a = a c i, suggesting the onset of a superfluid- 
like behaviour above this critical value of activity [llj . 
For a > a c i, the linearized approximation breaks down 
and the stress versus (average) strain rate curve ob- 
tained by numerical solution of the equations is nonlin- 
ear and nonmonotonic, as shown in Fig. [8] We empha- 
size that the flow profiles are always inhomogeneous with 
varying velocity gradients and director orientation. For 
a c \ < a < a C 2 the theoretical stress versus macroscopic 
(average) strain rate curve goes through the origin and 
exhibits a region of negative da/dj, that would in princi- 
ple be mechanically unstable. What would be measured 
in an experiment would, however, depend critically on 
details of the experimental procedure and the particular 
apparatus. To study the steady state rheology there are 
in general two natural classes of experiments: either (i) 
one tunes the stress a and measures the resulting strain 
rate 7 or (ii) one does a sweep through the values of strain 
rate 7 and measures the stress a. If the stress-strain rate 
curve is monotonic, the two procedures are expected to 
yield the same result. However, this is no longer the case 
as soon as the response exhibits nonmonotonicity. 

An important question, then, is what is the shape of 
the stress-strain rate curve that would be obtained exper- 
imentally for a > a c i in an experiment where one tunes 
the macroscopic strain rate 7. Several scenarios are pos- 
sible, as shown in Fig. [7] for a non-monotonic curve with 
maximum/minimum at ±cr m . 

(i) One scenario, suggested recently [TT] based on nu- 
merical studies in the proximity of the isotropic-nematic 
phase transition and for small value of the active stress 
a is the appearance of bulk shear bands accommodat- 
ing a range of macroscopic shear-rates at zero stress. 
This would correspond to the bulk stress-strain curve 
displayed in the top right frame of Fig. [7] and charac- 
terized as "superfluid" behavior. In the simplest picture 
the sheared suspension would separate in bands of con- 
stant and opposite strain rates, each with zero stress. 
For the systems studied here (deep in the ordered phase, 
either nematic or polar), we find that the equations of 
motion provide no mechanism for selecting a particular 
value of the stress plateau and are unable to find a stable 
stress-plateau at any value of \a\ < a m (including a = 0, 
see Fig. [7]). Furthermore we always find flow profiles 
with continuously varying gradients of fluid velocity for 
all values of macroscopic strain-rate 7 implying that the 
picture of two bands of constant strain rate would be at 
best an idealisation. 
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(ii) An alternative scenario that is observed in 
other driven systems, such as charge density waves in 
anisotropic metals [TTj and collections of motor pro- 
teins [IB], is shown in the bottom right frame of Fig. [7] 
In this case the system is expected to exhibit hysteresis, 
with regions that accommodate coexistence of a range of 
macroscopic strain rates, corresponding to the constant 
value ±<7o of applied stress. In general ao may coincide 
with o~ m or may be lower, with the system exhibiting 
"early swtching" . The width of the horizontal hysteretic 
region of the stress-strain curve decreases with increasing 
a. In this picture the particular steady-state behaviour 
observed will depend on the initial conditions and par- 
ticular flow history of each sample. 

(iii) Another possibility is that the system shows a 
yield-stress like behaviour with a yield stress ±a y whose 
sign is determined by the direction of the flow. The value 
of the yield stress could also be anywhere in the "unsta- 
ble" range of stress: a y < c m . 

(iv) Finally, there is one more possibility: that he 
theoretical curve would indeed be reproduced by an ex- 
periment which scanned through different values of the 
macroscopic strain rate. The theoretical curve has been 
calculated by fixing 7 and calculating the corresponding 
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FIG. 7: The top left frame display a typical theoretical stress- 
strain curve of a nematic active suspension in the region 
Qf c i < |q| < a C 2- The theoretical curve is obtained by tun- 
ing 7 and calculating the resulting a and exhibits a region 
of da/d A / < 0. The other three frames show three possi- 
ble experimental stress-strain curves obtained by tuning a 
and measuring 7 that could be consistent with the theoretical 
curve. The top right frame displays the "superfluid" scenario 
suggested in [TT|, with bulk shear bands accommodating dif- 
ferent macroscopic shear rates and zero net stress, so that 
the apparent viscosity of the system is simply zero. The bot- 
tom left frame shows a yield-stress like behaviour with a yield 
stress g v = a m . The last scenario is described in the bottom 
right frame and corresponds to a hysteretic stress-strain curve 
where the suspension can accommodate a range of macro- 
scopic strain rates maintaining a constant total stress ±<Trj. 
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FIG. 8: Stress-strain curves of a nematic suspension (fl = 
W = 0) obtained by numerical solution of the active hydro- 
dynamic equations for several values of a. a c \ = 0.219 and 
a C 2 — 0.877 for the parameters chosen in the numerical solu- 
tion. 




FIG. 9: Yield-stress <r c as a function of a for a nematic 
suspension (/? = w = 0) obtained by numerical solution of 
the active hydrodynamic equations. 



value of a under the assumption that there are variations 
in the director and flow field only in the gradient direc- 
tion (i.e. perpendicular to the plates). If this assumption 
is valid, every point on this curve does therefore represent 
a stable state corresponding to this procedure. 

For a > a C 2 the stress-strain curve intercepts the 7 = 
axis at a finite value a c — a(j = 0) of the strain rate. The 
active suspension has a nonzero spontaneous stress even 
in the absence of applied forces, as indeed observed in the 
spontaneous flow regime of an active suspension confined 
between two stationary no-slip planes. In other words, 
a finite force must be applied to the active suspension 
to keep it from sliding even at zero mean strain rate. 
This spontaneous stress a c is shown as a function of a in 
Fig. [9j The sign of the stress determines the direction of 
spontaneous flow. 

We now speculate on the possible behavior of the sys- 
tem for each of the scenarios sketched above as a goes 
through a C 2- The behavior is shown schematically in 
(i) In the superfluid scenario, the response of 
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the suspension to an applied macroscopic strain rate will 
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FIG. 10: Possible scenarios for the transition to the yield 
stress regime at a > a C 2- The non-monotonic curve obtained 
numerically is shown in the top left frame. In the superfluid 
scenario (top-right) the plateau at a — divides into two dis- 
connected branches terminating at a = ±<r c . In this case the 
yield stress is expected to grow monotonically from zero. In 
the yield stress scenario (bottom-left), there is already a non- 
zero stress at -y = and thus the yield stress simply continues 
increasing with no qualitative change in the behaviour at a C 2 ■ 
In the hysteretic scenario (bottom-right), the loop intersect 
the positive a axis at ±<to, with a c < ao < u m . 



show yield stress behavior. The system would smoothly 
go from the zero-stress plateau to a yield stress which 
increases from zero at a C 2- (ii) In the hysteretic scenario 
the minimum height of the hysteretic loop becomes 2a c 
i.e. a c < (To < <r m . (iii) In the yield-stress scenario the 
system already shows yield stress behaviour which con- 
tinues for a > a C 2- (iv) In the non-monotonic scenario, 
the non-monotonic stress-strain rate curve shows a jump 
at 7 whose magnitude increases from zero at a c i- 



V. DISCUSSION AND CONCLUSIONS 

We have studied the rhcological behavior of a thin film 
of polar and apolar active material. For weakly active 
systems, in the regime of the linear rheology, we have 
confirmed analytically the prediction of Hatwalne and 
collaborators [5] that activity can lower the linear bulk 
viscosity of tensile suspensions of swimmers as well as 
enhance the viscosity of contractile systems. We have 
shown that this result applies also for finite systems, in 
the presence of boundaries. 

An important new result of our work is the role of 
the shape of the active particles in controlling the The- 
ological behavior. We find a remarkable exact duality 
that holds in the regime where the stress-strain rate rela- 
tion is linear and shows that tensile (a < 0) rod-shaped 
flow-aligning particles (A > 1) are Theologically equiva- 
lent to contractile (a > 0) discotic flow-tumbling par- 



ticles (—1 < A < 0). This means that activity lowers 
the linear viscosity of both tensile, rod shaped particle 
and contractile, disc shaped particle suspensions, while 
it increases the linear viscosity of contractile, rod-shaped 
particle and tensile, discotic particle suspensions. 

For strongly active systems we find that the Theolog- 
ical response is intrinisically nonlinear. The regime of 
linear rheology at small strain rates vanishes beyond a 
critical value of activity. In this strongly active regime, 
we explore a number of possible scenarios for the non- 
linear rheology which include a "superfluid" phase with 
vanishing viscosity, hysteresis, yield-stress behavior and 
non-monotonic behavior. Our one-dimensional analysis 
does not, however, allow us to determine which of these 
scenarios is more likely. It is of course possible that allow- 
ing for variations of the director and flow field in higher 
dimensions or allowing for variations in the magnitude of 
the order parameter would yield a criterion for selecting 
one of the proposed scenarios. 
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Appendix A: Derivation of Eqs. (3) 

In this section we show some details of the derivation 
of the modified "passive" terms in the equation for the 
director field p in the polarized state, when fluctuations 
in the magnitude of the order parameters are neglected. 
The equation for the full vector order parameter P has 
the form 



[fit + v • V]Pi = XuijPj + Yhi + Vfi 



(Al) 



Eq. (All can be separated in two equations for the 
magnitude P = |P| of the polarization and its direction 
p = P/P, using 



d t p 



d t P = Pi d t Pi 
1 



(A2) 
(A3) 



where 6fj — Sij — P iPj is a transverse projection operator, 
with the result 



d t P = P(Xu ijPi p j )+T'f ]l +Th l] , 

[fi t + v • v]pi + ujijPj = xsfjUjkPk + nCr'/Z 



p 
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where we have defined 



ftii = P -h, h± = %h J , /n=p-f, /, ■ 

In the ordered state, fluctuations in the magnitude P of 
the polarization are overdamped and will be neglected. 
We can assume, on the other h and, to be deeply in the 
polarized state and that P = yf—a^/oi is constant. For 
simplicity we redefine the units so that P = 1. The 
condition P = const determines the longitudinal part /iy 
of the molecular field. This requires 

h w = -p[ r 7ii + ^uijPiPj] ■ 

The above expression can be now used to eliminate hn 
from the density j = T'h + T"f appearing at the right- 
hand side of Eq. (la). Expressing hi = Pih\\ + hf- and 
fi = Pif\\ + ft we obtain 



ji = Pi T"(l - - iXu klPk pm + T'ht + T"ft , (A4) 

where £ = (r') 2 /(Tr") is a dimensionless parameter and 
7' = r'/r. Similarly, the stress tensor becomes 



The longitudingal and transverse parts of the driving 
force fi are given by 



C -61 — B 3 
/ll = — 2 P • Vc - — 3 - p • V V • p , 

% c 



f ± =S T 



C 
r 2 



~3 d 3 c ^ V 'P 

Cn Cg 



Similarly, the transverse part of the molecular field is 
given by 



hi = Sjj 



Bi — B 3 
Co 



d jC + (K, - K 3 )d 3 V ■ p + K 3 X7 2 p 3 



Replacing the explicit expressions of h\\, hj~, /11 and ft 
in Eqs. (A4) and (A5), we finally obtain 



D{l-OPiPj-— Sfj 
co 



djC-j'XukiPkPiPi 



SijH- ^{Pihj 
1 

+ 2 



[pihf ~Pjht] - \piPjh\\ 



(A5) 



where w = r(Bj — B 3 ) is a velocity and we have neglected 
terms of second and higher order in the hydrodynamic 
fields. Finally, the reversible part of the stress tensor is 
given by 



ofj = -%H + XptPjPk 



Bi — B 3 
co 



d k c + {K x - K 3 )0 k V ■ p + K 3 V 2 p k 



Bi — B 3 
co 

Bi — B 3 
co 



(Pidjc + Pj dic) + {Kx - K 3 )( Pl dj + Pj di)V ■ p + Kaip^pj +p 3 V 2 p l ) 
{pidjC - pjdic) + (Ki - K 3 )(pidj - p 3 di)V ■ p + K 3 (p,V 2 p 3 - p 3 V 2 pi) 



\r'£piPj(DpkdkC + wpkdudipi) + —p l p 3 u kl p k p l 



Taking K\ = K 3 = K leads to the equations given in 
Sec. H 



[1] H. Gruler, U. Dewald, and M. Eberhardt, Eur. Phys. J. 

B 11, 187 (1999) 
[2] J. Toner and Y-h. Tu, Phys. Rev. Lett. 75, 4326 (1995); 

J. Toner, Y. Tu and S. Ramaswamy, Ann. Phys. 318, 

170 (2005). 

[3] R. A. Simha and S. Ramaswamy, Phys. Rev. Lett. 89, 
058101 (2002). 



[4] K. Kruse, J. F. Joanny, F. Jiilicher, J. Prost, and K. 

Sekimoto, Phys. Rev. Lett. 92, 078101 (2004). 
[5] T. B. Liverpool and M. C. Marchetti, Phys. Rev. Lett. 

90, 138102 (2003). 
[6] T. B. Liverpool and M. C. Marchetti, in Cell Motility, P. 

Lenz, ed. (Springer, New York, 2007). 
[7] R. Voituriez, J. F. Joanny and J. Prost, Europhys. Lett. 



10 



70, 118102 (2005). 
[8] L. Giomi, M. C. Marchetti and T. B. Liverpool, Phys. 

Rev. Lett. 101, 198101 (2008). 
[9] Y. Hatwalne, S. Ramaswamy, M. Rao and R. A. Simha, 

Phys. Rev. Lett. 92, 118191 (2004) 
[10] T. B. Liverpool, M. C. Marchetti, Phys. Rev. Lett. 97, 

268101 (2006) 

[11] M. E. Cates, S. M. Fielding, D. Marenduzzo, E. Orlandini 
and J. M. Yeomans, Phys. Rev. Lett. 101, 068102 (2008). 

[12] B. M. Haines, I. S. Aranson, L. Berlyand and D.A. 
Karpeev, Phys. Biol. 5, 046003 (2008). 

[13] A. Sokolov and I. S. Aranson, Phys. Rev. Lett. 103, 
148101 (2009). 

[14] B. M. Haines, A. Sokolov, I. S. Aranson, L. Berlyand and 



D.A. Karpeev, Phys. Rev. E 80, 041922 (2009), 

[15] A. Baskaran and M. C. Marchetti, Proc. Natl. Acad. Sci. 
USA 106, 15567 (2009). 

[16] The parameters Pi and a have different dimensions from 
those used in [6] as here we have incorporated an addi- 
tional factor of f 2 in their definition. 

[17] A. Maeda, M. Notomi and K. Uchinokura, Phys. Rev. B 
42, 3290 (1990). 

[18] F. Julicher and J. Prost, Phys. Rev. Lett. 75, 2618 
(1995). 

[19] We neglect here convective nonlinearities in the Navier- 
Stokes equations that are unimportant on the long time 
scales of interest. 



